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$H ' We apply the time-dependent variational principle, the nuclear field the- 

ory, and the boson expansion method to the Lipkin model to discuss anhar- 
monicities of collective vibrational excitations. It is shown that all of these 



' approaches lead to the same anharmonicity to leading order in the number 

> ' 



of particles. Comparison with the exact solution of the Lipkin model shows 
that these theories reproduce it quite well. 



^ , I. INTRODUCTION 



The existence of multiphonon states in the nuclear spectrum of excitation has been 



^ i predicted since the introduction of collective models Examples of low- lying nuclear 
vibrational states have been known for many years in nuclear spectra and are still being 
actively investigated with new generation of detectors; in particular, two-phonon multiplets 
and some three-phonon states based on the low-lying collective quadrupole and octupole 
modes have been found Recently, it has also been beautifully demonstrated that 

multi-phonon excitations of these low-lying collective vibrations strongly influence heavy- 
ion fusion reactions at energies near and below the Coulomb barrier [^], through the so- 
called fusion barrier distribution analysis pO|JTl[] . It was pointed out that anharmonicities 



of vibrational excitations can significantly alter the shape of fusion barrier distribution and 
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that thus sub-barrier fusion reactions offer an alternative method for extracting the static 
quadrupole moments of phonon states in spherical nuclei |T^|T^. 

In the past 15 years, evidence has been collected for two-phonon giant resonances as 



well [14|. This evidence stems from heavy-ion reactions at intermediate energy |15,16 



pion-induced double charge exchange reactions |T^, and relativistic heavy-ion reactions via 



electromagnetic excitations, in particular the excitation of the double giant dipole resonance 
(DGDR) [|l^. In the experiments of the last type, a puzzling problem has been reported 
|T9|^2|. Although the experimental data show that the centroid of the DGDR is about 
twice the energy of the single phonon resonance, theoretical calculations which assume the 
harmonic oscillator for giant resonances considerably underestimate the cross sections for 
double phonon states. In connection with this problem, anharmonic properties of giant 
resonances are attracting much interests p3|^9[ . 

Recently Bertsch and Feldmeier applied the time-dependent variational principle 
(TDVP) [^] to large amplitude collective motions and discussed anharmonicities of various 
giant resonances |^. One of the advantages of their approach is that solving the resultant 
equations and estimating the degree of anharmonicity are quite simple. They found that the 
relative frequency of the double phonon state scales as /S.uj/uj ~ A~^/^^ A being the atomic 
number of a nucleus. 

Earlier, Ponomarev et al. [0 noted that this quantity scales as y4~^/^ in the nuclear 
field theory (NFT) [^], the same as the NFT result for the octupole mode Reference 
[|l|] also remarks that a liquid-drop estimate gives an A~^/^ dependences, implying that 
quantal effects are responsible for the difference. Both NFT and TDVP are quantal theories 
giving different results in refs. and p3[, so differences must be due either to inadeqate 



approximations or differences in the underlying Hamiltonians. We therefore undertook to try 
both methods on a solvable Hamiltonian. This will test the reliability of both methods, and 
if both give correct results, the disagreement is very likely attributable to the Hamiltonian 
assumptions. 

The time-dependent variational approach was recently applied also to dipole plasmon 



resonances of metal clusters p3[. In Ref. [|33|, it was shown that the anharmonicity of 
dipole plasmon resonances is very small and scales as Acu/cj ~ A~^^'^, which is identical 
to the result of the time-dependent variational approach for nuclear giant resonances. On 
marked contrast, Catara et al. claimed a large anharmonicity of plasmon resonances using a 
boson expansion method , which contradicts both to the result of Ref. [^] and to recent 
experimental observation This fact also motivated us to compare these two methods, 
together with the nuclear field theory, on a solvable model in order to clarify the origin of 
the discrepancy. 

As an aside, we mention also that there is a size-dependent anharmonicity in quantum 
electrodynamics, considering the photon spectrum in a small cavity. In QED the only dimen- 
sional quantity is the electron mass mg, and the photon-photon interaction is proportional 
to Thus on dimensional grounds the relative shift of a two photon state in a cavity of 
side L scales as Auj/uj ~ l/m^L'^. Considering that sizes of a system scale as i? ~ A^/^, the 
results of refs. [p3| , p3| are /S.uj/uj ~ identical to QED. 

In this work we compare the nuclear field theory, the time-dependent variation approach, 
and the boson expansion method on the second Lipkin model, eq. (2.5) in ref. The 
model is finite-dimensional and can be solved exactly. It has widely been used in literature 
to test a number of many-body theories |]37|- ^] . 

The paper is organized as follows. In Sec. II, we first solve the model in the random phase 
approximation (RPA) and discuss the harmonic limit. In Sect. Ill we derive the collective 
Hamiltonian using the TDVP. We requantize the collective Hamiltonian and discuss the 
deviation from the harmonic limit. Numerical calculations are performed and compared 
with the exact solutions. In Sec. IV, we use the nuclear field theory as an alternative. There 
we see that it gives the same result as the TDVP, to leading order in the dependence on 
the number of particles. In Sec. V, we compare those results with the boson expansion 
approach. It will be shown that it leads to the identical result to the TDVP and the NFT. 
Finally we summarise the paper in Sec. VI. 
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II. HARMONIC LIMIT 



1=1 

N 



Lipkin, Meshkov and Glick |36] proposed two Hamiltonian models to describe N particles, 
each of which can be in two states, making 2^ states in the basis. Using Pauli matrix 
representation for the operators in the two-state space, the second model has the Hamiltonian 

1 V 

H = -ecr^ - —cr^^aa:. (1) 

The first term is the single-particle Hamiltonian with an excitation energy e, and the second 
term is a two-body interaction. The quasi-spin operators and ax are given by 

N 

^{(^uO'ii — %i(^oi) , (2) 

i=l 

N 

= '^{(Ai^^oi + oliaii), (3) 

i=l 

respectively. [an) and ajj (aoj) are the creation (annihilation) operators of the i-th 
particle for the upper and the lower levels, respectively. For small V, the Hartree ground 
state |0 > is the fully spin-polarised state with matrix elements given by 

< 0|a,|0 >= -NS,,,. (4) 

A suitable basis for the exact diagonalization of H is the set of eigenvectors of the angu- 
lar momentum operators and Jz with J = N/2. Then the dimension of the matrix 
diagonalization is reduced from 2^ to -|- 1 , making the numerical problem very easy. 

Before going to the anharmonicity, we note that the harmonic limit is obtained by solving 
the RPA equations. This was carried out in ref. for the first Lipkin model. The RPA 
frequency for the second Lipkin model, Eq. (p, is obtained in the same manner. Setting 
the RPA excitation operator as 

N \ / N 



(5) 



the RPA equations read 
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A, 



(x\ 



v° -1/ 



(6) 



where A and B are given by 



1 



^=4^<0|[a_,[iJ,a+]]|0>=e(l-x), 



B 



4iV 



<0|[a_,[/f,a_]]|0>=ex, 



respectively. We have defined o"_, cr+, and x as 



1=1 



1=1 



X = l^(iV-l)/e, 



(7) 
(8) 

(9) 

(10) 
(11) 



respectively. From these equations, the RPA frequency and the amplitude of the forward 
and the backward scatterings are found to be 



X 
Y 



{A-rB){A-B) = e^\-2x. 
uj + e 



UJ — e 



(12) 
(13) 

(14) 



respectively. 

Fig. 1 compares the exact solution for the excitation energy of the first excited state 
with the RPA frequency given by eq. (p!2D. As a typical example, we set the strength of 
the interaction V/e to be 0.01. The solid line shows the exact solutions for this particular 
choice of the Hamiltonian parameters, while the dashed line shows the RPA frequency. 
AtA^ = 51, x = 0.5, the RPA frequency becomes zero and the system undergoes phase 
transition from spherical to deformed. 

Figure 2 is the same as Fig. 1, but for a fixed x- We set x to be 0.25, which corresponds to 
the isoscalar giant quadrupole resonance. We find significant deviation of the RPA frequency 
from the exact solution for small values of A^, suggesting large anharmonicities. We discuss 
now the deviation from the harmonic limit. 



III. TIME-DEPENDENT VARIATIONAL APPROACH 



The time-dependent variational approach has been apphed to the first Lipkin model by 



Kan et al. in ref. but has never been applied to our knowledge to the second model. In 
keeping with the procedure of ref. we postulate a time-dependent wave function of the 
form 

\a(5) = exp(iQ;(t)(T2;) exp(i/3(t)crj^)|0). (15) 



The motivation for this ansatz appears in ref. [^. The operator in the first term is the one 



that we wish to evaluate in the transition matrix elements. The operator in the second term 
is obtained by the commutator with the Hamiltonian, 

[H,aJ^ = itay. (16) 

The Lagrangian is given by 

C = -a{f3\a.M - (17) 

We reduce this with the help of the identity 

e-^'^^Vj-e*'"'^ = cos 26 aj + sin 29 a ■ (I x j), (18) 

where i ^ j are Cartesian indices of the Pauli matrices. For example, the bracket in the 
first term of eq. (0) is reduced as 

{f3\a,\/3) = (0|e-^'^-^a,e*'^'''5|0) =cos2/?(0|a,|0) -sin2/5(0|a,|0). (19) 

The first term in the Lagrangian is 

-a{P\a^\p) = -Nasm2(3. (20) 

The second term is 

- {ap\H\ap) = eyCos2a cos 2/? + V'y (cos^ 2/? + A^sin^ 2/5 ). (21) 



The Lagrangian may then be expressed as 

C = -Na sm2p + e^ cos 2a cos 2/3 + V"^ (cos^ 2P + N sin^ 2/3 ) . (22) 

The first Lagrangian equation is d/dtdC/da — dC/da = 0. It reduces to 

/3 = ^ sin 2a. (23) 

Similarly from the second Lagrange equation, d/dt dC/d(3 — dC/df3 = we obtain 

a cos 2/? + I cos 2a sin 2(] - {N - 1)V cos 2/3 sin 2/3 = 0. (24) 

Next let us linearize and see what the RPA frequencies would be. The linearized equations 
are 

p = ea (25) 

a + {e~2V{N -l))f3 = 0. 
The equation for the frequency reads 

uj"^ = e"^ -2{N -l)Ve = e\l -2x) (26) 



in agreement with the result of eq. (p!2|). 

A Hamiltonian corresponding to our Lagrangian can be seen by inspection, comparing 



eq. (p2D to the form 

C = qp-n{p,q). (27) 

Equation (p2|) is already of this form with e.g., p = sin 2/3, q = 2a. The Hamiltonian is 
then given by 

/ : — ^-7 ..A^ 



nip, q) = -ey cosg^/l - (2p/Ar)2 - V^-[{1 - {2p/Nf) + N{2p/Nf). (28) 

We now expand Ti in powers of q and p up to fourth order. Dropping the constant term, 
the expansion has the form 



with coefficients 



U=^ + ^q^ + aq^ + bp' + cqV (29) 
2m 2 



14 (.-2m- 4 (30) 
k = e-, (31) 



^=4' (33) 
c = (34) 



48 
e 

Yn' 

Note that we recover the hnear frequency, eq. (|T^, immediately from uj^ = k/m. 

In ref . , the anharmonicity was determined by requantizing the Hamiltonian with the 
Bohr-Sommerfeld condition, 



/ pdq = niT, (35) 

JqO 



where p and q satisfy Ti.{p, q) = E and go and qi are the endpoints of the motion at energy 
E. However, here we find it more convenient to use the equivalent formula 

/ dpdq = 2nTC. (36) 

In the same sense as the expansion of the Hamiltonian Ti.{p,q) as done in eq. (p9D, we 
apply eq. ( |36D iteratively. We ffist consider only the harmonic part of the Hamiltonian 
and transform the integration region to a circle. We also use polar coordinates and write 
p' = p/ y/m = r sin 6,xq = \/kq = cos 9. The radius of the circle is then Tq = \/2E and the 
harmonic approximation gives 

r Ira f 1 f^o E 

/ dpdq = J— / dp'dq' =- dO rdr = 2-K-. (37) 
Jh<e V k Jh<e u Jo Jo uj 

The nonlinearity can now be treated as a perturbation. To lowest order in the quartic terms, 
the radius to the boundary surface is given by 
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rl ^ rl - 2rl ( cos^ Q + m% sin^ ^ + ^ sin^ Q cos' Q \ . (38) 

k J 

This integral is easily evaluated with the result 

— / dvdq^ ---[-] 3— + 36m' + ^ • (39) 



27r 77Y<£; c<j 2 \uj J \ k"^ k 

Inserting the parameters from eqs. (|2P|)-(R¥), the anharmonic term can be expressed 



UJ /E^ ' 



8eA^ \uj 



) (-l + 3(eM^-2(e/.;)'). (40) 



Note that if there is no interaction, uj = e and the anharmonicity vanishes. This is rather 
remarkable; the Hamiltonian in this case is the first term in eq. (pS]), which looks nonlinear. 
But the solution of the equations of motion are independent on excitation energy. It is 
not a harmonic oscillator spectrum, however, because the energy is bounded. These two 
properties correspond exactly to the quantum spectrum of the operator eJz- 
We next quantize the above action to get 

8eN' 



En = nuj + (-1 + 3(e/cu)4 - 2{e/ujY) . (41) 



Taking the second difference, this yields an anharmonicity of 

,2 , 0,.^3 



The exact value of anharmonicity A^'^^E is compared with the value obtained from 



eq. (^21) in Fig. 3. We can see that the time-dependent variational principle works very 
well. 



IV. NUCLEAR FIELD THEORY (NET) APPROACH 

The NFT is a formulation of many-body perturbation theory with vibrational modes 
summed to all orders in RPA. Its building blocks are RPA phonons and the single particle 
degrees of freedom which are described in the Hartree-Fock approximation. The coupling 
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between them is treated diagramatically in the perturbation theory. For the Hamiltonian 
(1^), the effective NFT Hamiltonian is given, to the lowest order, by 

Hnft = ^ea, + uO^O + Hp^. (43) 

The first term in the Hj^ft describes single-particle spectrum. In writing down this term, 
we have used the fact that, for small value of V, the excitation energy in the HP is given 
by e and that the creation and the annihilation operators for the HF levels are the same as 
those for unperturbative levels. The second term describes the RPA phonons, with u and 
given by eqs. (|12|) and (^), respectively. The particle- vibration interaction Hp^ in eq. 
(1^) is given as 

Hp, = ~A{0^ + 0)a,, (44) 
where the coupling constant A is given by 



x' being NV/e. This Hamiltonian is constructed by replacing in the two-body interaction 
in the original Hamiltonian (|I|) as in Ref. In general, there is also a residual interaction 



among particles and holes, but it does not contribute for the Hamiltonian ([l|) to the lowest 
order. 

Each graph in the NFT contributes to a given order in but to all orders in x'- Since 
the microscopic origin of the RPA phonon is a coherent sum of particle-hole excitations, 
bubble diagrams have to be excluded when one calculates physical quantities in the NFT 
||32|| . To the zero-order (in 1/A^), the phonon energy coincides with that in the RPA given by 
eq. (|12|). The anharmonicity begins with the leading diagrams, which are shown in Fig. 
4. These diagrams are called "butterfly" graphs (see also refs. |l|, p2| , ^ ). For each diagram 
shown in fig.4, there are 5 other diagrams which are obtained by changing the direction of 
the phonon lines. As already said, for the Hamiltonian (^ there is no diagram of order 
involving residual interaction among fermions. 
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The contribution from each diagram is most easily evaluated by using the Rayleigh- 
Schrodinger energy denominator, which are more suitable in the lowest order expansion 
The four graphs in Fig. 4 have identical contributions, each given by 



In this equation, the minus sign appears because of the crossing of two fermion lines 
By summing up the contributions from all diagrams, we obtain 

To compare with eq. (^), note that x' = X the leading order of With this substi- 
tution, the two results are identical. 



V. BOSON EXPANSION APPROACH 

In the boson expansion method, each fermionic operators is replaced by corresponding 
operators which are written in terms of boson operators. There are several prescriptions 



to carry out the mapping from the fermionic to the bosonic spaces |^9|. Here we follow 



Refs. p9| , p| , ^ which discussed the anharmonicities of the Lipkin model using the Holstein- 
Primakoff mapping. In this mapping, fermionic operators are mapped to bosonic operators 
so that the commutation relations among operators are preserved. The quasi-spin operators 



in the present two-level problem are then mapped as 



hi 



a^-.2VNB^Jl-^ (48) 



a_ ^ 2 ViVy 1 - ^j^B (49) 
a,^-N + 2B^B, (50) 

where the operators B and 5^ satisfy the boson commutation relation, i.e., [5,5^] = 1. 
The Hamiltonian in the boson space which corresponds to the Lipkin Hamiltonian (|I]) is 
therefore obtained as 
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+VB^B + ^ (fiT^t + + VB^B^BB + ^ (b^ B^ B^ B + B^BBB) , (51) 
to the first order in 

A truncation of tlie expansion up to tlie leading order of the corresponds to the 
RPA which we discussed in Sec. II. To this order, the boson Hamiltonian (|5T|) is given by 



H"^^' = e[l- —^j B^B - [B^B^ + BB) . (52) 

As is well known, this Hamiltonian can be diagonalised by introducing a transformation 

= XqO^ + YoO (53) 

B = XoO + YoO\ (54) 

and imposing 

[hS\0^]=ujO\ (55) 
with a condition Xq —Yq=1. The frequency uj then reads 



uj = e^l- 2x', (56) 
x' being NV/e, which was introduced in the previous section, together with 

(5T) 

r^^^^±^, (58, 
x„y„ = g. (59) 

The frequency uj coincides with that in the RPA given by eq. (p!2| ) to the leading order of 
1/N. 

As in the nuclear field theory approach, the anharmonicity begins with the next order of 
the expansion. In terms of RPA phonon creation and annihilation operators defined in 



eqs. ( |53D and (|5^) , the boson Hamiltonian (|51|) can be rewritten as 
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+if4o (O^^ + O") + Hu [O^^O + 0^0^) + H220^^0\ (60) 

where the first term is the leading order of the expansion given by eq. (^) and the 
rest are the higer order corrections. The coefficient -^225 for example, is given by 

i^22 = V{X^^ + SX^Fo + 4X2^^2 + 3XoF„3 + Y^) = ^ - . (61) 

In order to estimate the degree of anharmonicity, we use the perturbation theory. The first 
order perturbation gives the energy of the one and the two phonon states of 

Ei = u + Hn, (62) 
E2 = 2{uj + Hu) + 2H22, (63) 

respectively. Taking the second difference, the anharmonicity reads 

A<=)£ = 2ft, = |g (1 - I) , (64) 

which is identical to the result of the variational approach given by eq. (|4^) as well as that 
of the nuclear field theory, eq. (^7[). 



VI. CONCLUSION 

We have shown that the nuclear field theory, the time-dependent variational principle, 
and the boson expansion method give identical leading-order anharmonicities for the Lipkin 
model, and that the formulas agree well with the exact numerical solution. The anhar- 
monicity is inversely proportional to the number of particles in the system, when the other 
parameters are fixed to keep the harmonic frequency the same. This clarifies the origin 



of the conflicting results for the A-dependence of the anharmonicity obtained in ref. |23 



and [32|. In ref. 123] the time-dependent method was applied to a Skyrme-like Hamiltonian 



involving all A nucleons, and the result was A(2)E oc fiuj)/A. In ref. fS^, the Hamihonian 
was restricted to a space of a single major shell for particle orbitals and similarly for the hole 
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orbitals. Since the number of particles in the valence shell increases as A^/^, the result was 
A^'^'^E oc /(co')/y4^/^. Finally, it should perhaps be emphasized that both methods predict 
that the anharmonicity is very small for giant resonances: both A^^'^ and A are large num- 
bers. This need not be the case for low-lying collective vibrations. The NFT can produce 
large effects when there are small energy denominators. Low-lying excitations are difficult 



to describe with a simple ansatz like eq. (|T5D, so the time-dependent variational principle is 
not easily applied. Clearly this is an area that should be explored further. 

As for the discrepancy between the time-dependent variational approach and the boson 
expansion method concerning anharmonicities of plasmon resonances of metal clusters, our 
study showed that the origin of the discrepancy should not be ascribed to the method used 
to solve the problem. The origin of the discrepancy, therefore, is not traceable at moment 
and further studies may be necessary in order to reconcile the discrepancy. 
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Figure Captions 



Fig. 1: Energy of the one phonon state as a function of the number of particle. The sohd 
hne is the exact solution of the Lipkin model, whilist the dashed line is obtained in the RPA. 
The strength of the interaction Vje is set to be 0.01. 

Fig. 2: Same as fig. 1, but for a fixed x parameter, which is set to be 0.25. 
Fig. 3: Anharmonicity of the double phonon state as a function of the number of particle 
for X = 0.25. The dashed line is an approximate solution given by eq. (]4^), while the solid 
line is obtained by numerically diagonalising the Hamiltonian. 

Fig. 4: The diagrams of order 1/A^ contributing to the anharmonicity of the double phonon 
state for the Hamiltonian (|l|). The wavy lines represent phonon propagations, while the 
particles and the holes are depicted by the allowed lines. 
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